Threshold detachment of negative ions by electron impact 
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The description of threshold fragmentation under long range repulsive forces is presented. The 
dominant energy dependence near threshold is isolated by decomposing the cross section into a 
product of a back ground part and a barrier penetration probability resulting from the repulsive 
Coulomb interaction. This tunneling probability contains the dominant energy variation and it 
can be calculated analytically based on the same principles as Wannier's description for threshold 
ionization under attractive forces. Good agreement is found with the available experimental cross 
sections on detachment by electron impact from D~ , 0~ and B~ . 



PACS numbers: 3.65Sq, 34.80D, 34.10+x 

Storage ring based experiments on threshold detach- 
ment from the deuteron (D~) and the oxygen (0~) neg- 
ative ions by electron impact |l],|2|,||], and recently also 
from B~ 0, have stimulated the theoretical interest in 
the mechanism and the quantitative description of this 
process |||6|,0,f|[§ . It is a fundamental question how 
threshold detachment proceeds since for very low energies 
the impacting electron does not even reach the atom be- 
cause it is repelled by the loosely bound electron. Early 
theoretical work on this problem tried to describe the 
process by asymptotic properties of the wavefunction for 
the two electrons in the continuum after the collision jl0| , 
following the spirit of Wigner's treatment for two-body 
break up (TT). However, the predicted cross section agrees 
purely with the experimental results. 

Some recent theoretical treatments, following another 
idea of the early days (l2| , emphasize the importance of 
tunneling contributions, either by treating the impacting 
electron as a constant perturbing electric field Q, or by 
merging a quantum and a classical description As- 
tonishingly good agreement with the experiment, even at 
low energies near threshold, comes from a coupled chan- 
nel calculation in the impact-parameter formalism where 
a classical trajectory is used for the relative motion of 
target and projectile electron and the electron to be de- 
tached is described quantum mechanically . These re- 
sults, at least the shape of the cross section, depend little 
on the polarization potential used, as Lin etal. emphasize 
. Results of similar accuracy have have been reported 
using a lowest order distorted-wave scheme, however, in 
contrast to , with a sensitive dependence on the polar- 
ization potential ||. 

Without a full calculation of all electrons, one cannot 
avoid to use parameters in one or another way, either 
directly in the simpler models Jl|,||, or indirectly in the 
more involved calculations modeling polarization poten- 
tials for the loosely bound electron 



The theoretical work so far remains inconclusive con- 
cerning a dominant mechanism of threshold detachment, 
and the reason for the seemingly contradicting findings 
concerning the robustness of the results with respect to 
changes in the polarization potential is unknown. 

A successful description of near threshold detachment 
focusing on threshold properties should naturally depend 
very little on details of the polarization since the long 
range repulsion between target and projectile electron 
dominates. Moreover, such an approach should uncover 
a mechanism for threshold detachment and thereby clar- 
ify the issue of robustness with respect to different polar- 
ization potentials. 

In the following we will show that threshold detach- 
ment by electrons can be described with the same tech- 
nique which has lead to the successful (and purely classi- 
cal) description of threshold ionization under long range 
attractive Coulomb forces, pioneered by Wannicr |13|. 

However, in order to learn how to deal with repul- 
sive Coulomb forces, one must go back to a semiclassical 
formulation of threshold ionization and analyze the rea- 
son why Wannier's classical treatment was appropriate. 
Semiclassically, one may write the scattering amplitude 
in the form fll41 



(1) 



where the sum runs over all scattering orbits j which 
contribute with the weight y/Vj. The phase contains the 
Maslov indices Vj III]] and the action $>j along the orbit 
which may be expressed as 

<P J (E)=<j )j (E)E-^ (2) 



where <f>j(E — > 0) = const. |Tj]. This special form is a 
consequence of the homogeneous Coulomb interaction. It 
is crucial for the justification of the classical treatment 



1 



since E — ► 0, i.e. approaching threshold, and Ti — > have 
the same effect in Eq. ([!]). If $j is real, which is the 
case for all classically allowed trajectories, one arrives by 
stationary phase approximation (for E — > or h — » 0) at 
the result 



(3) 



which sums all individual contributions Vj of the trajec- 
tories to the classical cross section ocl- 

Looking for the dominant energy dependence of a(E — > 
0) we decompose the cross section into 



a(E) = <j b (E)P(E) 



(4) 



where as(E) is a smooth back ground cross section with 
<jb(E — > 0) =const. Wannier showed that the domi- 
nant energy dependence P{E) — P* is contained in a 
single fixed point orbit J = * | [L3| . Formally, this or- 
bit represents an outgoing trajectory with fixed angle 
6* = 7r between the two electrons and symmetric dis- 
tances fx = ^2 of electron 1 and 2 from the core. It 
is convenient to use hyperspherical coordinates with an 



overall radius r 2 



of the system and the hyper- 



angle defined by tana = ri/f2- The symmetric escape 
orbit corresponds to increasing hyperradius r and fixed 
a* = 7r/4. The potential energy of the two electrons in- 
teracting with a core of charge Z can be written in the 
form of a Coulomb potential with an angular dependent 
charge, V = C(a, 0)/r. For Z > 1/4 the potential at the 
fixed point is with C* = C(a*,6*) < attractive. Hence 
the relevant threshold orbit is classically allowed with a 
real action <!>» and, as sketched above, the semiclassical 
scattering amplitude leads for E — > to the classical 
cross section with dominant energy variation of the form 



Pcl(E) = (E/E ) 



(5) 



as derived by Wannier [^3[ . 

On the other hand, for a fixed point charge C» > 
the Coulomb interaction is repulsive. Then, the relevant 
threshold orbit is classically forbidden and represents a 
tunneling trajectory with imaginary action $* = ir*. In 
this case the semiclassical cross section does not reduce 
to the classical one in the limit E — ► 0. Rather, its major 
energy dependence results from a tunneling mechanism 
which produces a Gamow factor 



P(E) = exp[-2T*(E)/ti\. 



(6) 



Clearly, the threshold cross section is through Eq. (|) h- 
dependent. Nevertheless, the important dynamical quan- 
tities, namely the tunneling action T*, is still given clas- 
sically. 

In the traditional description of classically allowed 
threshold fragmentation of charged particles the initial 
configuration is unimportant - the energy dependence of 
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FIG. 1. Sketch of tunneling threshold dynamics on the 
fixed point manifold with potential C*(r) from Eq. (^). The 
classically allowed incoming and outgoing trajectory on the 
respective energies Ei and Ef is shown (dashed lines) as well 
as the tunneling part (solid thick line) which determines the 
threshold fragmentation probability. 

the cross section is completely determined by the stability 
of the final fragment configuration. This stability enters 
P* of the escape orbit. That only the fixed point orbit is 
relevant close to threshold is justified by the fact that all 
available energy (which approaches zero for E — > and 
r — * oo) must be put into the radial degree of freedom 
r in order to fragment the system. Hence, the system 
evolves asymptotically in a frozen configuration where 
neither its geometrical shape (8 = 9*) 1 nor the relative 
interparticle distances ri/ra = tana* change. Moreover, 
due to the Coulomb scaling properties, any partial wave 
with angular momentum L reduces in scaled coordinates 
to an S'-wave since the scaled angular momentum reads 
L = L\fE p4| . Therefore, only the S'-wave has to be con- 
sidered which remains also valid in the case of a repulsive 
Coulomb force. Finally, for two escaping electrons, the 
fixed point configuration 9* — 7r and a* — n/4 remains 
the same for all charges of the core including the limit 
Z = which applies to the neutral atom for our prob- 
lem of electron detachment. Hence, the radial motion on 
the fixed point manifold is governed by the Hamiltonian 
(atomic units are used unless otherwise stated) 



H* 



2 r 



(7) 



where the effective charge C* = 2 x l 2 results from the 
evaluation of the electron-electron repulsion 1^=17^1 — 



''2 



at the fixed point. 



For each energy E = H* we can calculate the tun- 
neling action T*(E) entering Eq. (H) from the imaginary 
momentum p = (— P 2 ) 1 / 2 of Eq. (R), 



r* = 



pdr. 



(8) 
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FIG. 2. Detachment cross section by electron impact as a 
function of excess energy for B~ (circles) from [Q, 0~ (tri- 
angles) and D~~ (diamonds) from 111. The solid lines are the 
cross sections from Eq. (^J). 

The integration limits are the outer turning point r t 
where the orbit becomes classically allowed, p(r t ) = 0, 
and a starting point rj, see Fig. 1. In contrast to thresh- 
old fragmentation under attractive Coulomb forces tun- 
neling threshold fragmentation depends on the initial 
configuration, at least as far as the value of in Eq. (|^) 
is concerned which will influence shape and magnitude of 
P{E) in Eq. (§). 

In a very crude approximation one could put ri = ar- 
guing that the the electronic momentum transfer requires 
the recoil to be absorbed by the nucleus and its position 
is where the outgoing electrons should start. However, 
in the light of the (small) tunneling probability which 
determines threshold detachment according to Eq. @ 
close to E — this is certainly too crude. For small ex- 
cess energy the projectile electron impacts roughly with 
the binding energy / which is of the order of 1 eV. Re- 
pelled by the loosely bound electron the projectile will 
never reach « at this low impact energy. More real- 
istically, one can approximate ri by the classical turning 
point of the incoming electron, projected onto the fixed 
point manifold whose dynamics is specified by Eq. ([?]). 
Hence, to determine this turning point of the incoming 
electron we put P r — in Eq. (0) at the incoming electron 
energy of J5» = E + 1 to yield 



n = c./{E + i). 



(9) 



The initial momentum of the outgoing electron pair 
p(r;) = \/2I follows from the Hamiltonian Eq. @) on the 
final energy surface Ef = E. The situation is sketched 
in Fig. 1. Using Eq. (||) and Eq. @ the threshold de- 
tachment probability Eq. (|^) reads in dimensionless units 
explicitly 



FIG. 3. Experimental detachment probabilities, obtained 
by dividing the cross section by gb{E) from Eq. (p"l|). The 
coding of the data is as in Fig. 3. Theoretical P(E) from Eq. 

<§>• 



P(E) = exp 
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(10) 



where a = 1/137 is the fine structure constant, m e c 2 = 
511keV is the rest mass of the electron, and C — C* is the 
repelling charge of the two electrons on the fixed point 
manifold in units of e, see Eq. (Q). One can cast Eq. jl0| ) 
into a more familiar form of atomic units by noting that 
m e c 2 /a 2 = e 2 /ao = 27.21 16eV is just the atomic energy 
unit. Clearly, the tunneling mechanism breaks the scal- 
ing invariance of P(E) for different systems characterized 
by different ionization potentials I since P(E) does not 
only depend on E/I but also on m^c 2 /E. This is one 
of the major differences compared to Wannier's classical 
result (Eq. (||)) for threshold ionization under attractive 
Coulomb forces. 

Different P(E) are shown in Fig. 3 with solid lines cor- 
responding to detachment from the ions B~ , D~ , and 0~ 
respectively. The "experimental" tunneling probabilities 
are extracted by fitting the experimental cross sections 
(Fig. 2) to Eq. (|) with 



a B (E) = a /(b + E/I), 



(11) 



where 00,60 are fitting parameters. The o~b(E) obtained 
in this way are shown in Fig. 4 for completeness and 
exhibit the expected monotonically decreasing behavior. 

As a final support for the analytical P(E) from Eq. 
(|io| ) we have fitted the experimental cross sections with 
(To, bo and / as free parameters. The result for / was 
0.297 ± 0.008,0.79 ± 0.03 and 1.58 ± 0.04eV which is 
close to the accurate values of 0.28,0.75 and 1.46eV for 
B~,D~, and 0~, respectively. 
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while Pcl(E) (Eq. (|J)) under attractive Coulomb forces 
is scale invariant. Yet, P{E) for threshold detachment 
can be described semiclassically due to the dominant (re- 
pulsive) Coulomb interaction which ensures through its 
scaling properties that E — > also means Ti — > (see 
Eq. (^)). The same scaling properties also reduce the 
dominant energy dependence of all partial waves to that 
of L = 0. Therefore, P{E) can be determined from the 
S-wave only, as it has been done in the present work. 

I would like to thank L. H. Andersen for making avail- 
able the experimental data electronically, in the case of 
B~ even prior to publication. Financial support from the 
DFG under the Gerhard Hess-Programm is gratefully ac- 
knowledged. 



FIG. 4. Back ground cross sections ob{E) for B (solid), 
D~ (dashed) and 0~ (dotted), see text. 

The present description differs from various published 
tunneling models approximating in one or another way 
the actual electron motion by tunneling. In the present 
treatment, only the dominant energy dependence of the 
cross section is derived from a fixed point orbit which 
represents a tunneling trajectory. However, this trajec- 
tory does not correspond to a true, physical two electron 
orbit. Rather, it is a stationary point solution for h — > 0, 
in complete analogy to Wannier's solution for the clas- 
sically allowed case of attractive forces. This stationary 
point calculated in the limit E — > does not depend at all 
on the polarization potential. Only the binding energy of 
the target electron enters P{E) through n as defined in 
Eq. (^) from the turning point of the incoming trajectory. 
It is exactly this element which is similarly contained in 
the calculation of Ref . [Q. Hence, this impact parameter 
calculation captures an essential feature of the threshold 
detachment dynamics making the whole calculation ro- 
bust against details of the polarization potential. These 
details will influence on the other hand the background 
cross section <7b{E) much more strongly. The distorted 
wave calculation |8| by nature approximates the thresh- 
old region from an expansion of the high energy limit 
which is much more sensitive on details of the (shorter 
range) polarization potential. 

In summary, separating the rapidly changing detach- 
ment probability P(E) from the background cross section 
ctb(E) we have shown that threshold fragmentation un- 
der asymptotic repulsive Coulomb forces can be treated 
on the same footing as the well established threshold ion- 
ization under attractive Coulomb forces. In contrast to 
the classical result for attractive forces, threshold detach- 
ment of negative ions by electrons can be interpreted to 
proceed via quantum mechanical tunneling of the out- 
going electron pair. This implies a breaking of the scale 
invariance of P{E) with respect to energy since P(E) de- 
pends on m e c 2 /E irrespectively of the target properties 
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